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Q ; ABSTRACT 

. We study the possibility of generating magnetic fields during the evolution of elec- 

' ' tron, proton, and photon plasma in the pre-recombination era. We show that a small 

, magnetic field can be generated in the second order of perturbation theory for scalar 

^ ' modes with adiabatic initial conditions. The amplitude of the field is ^ lO"'^" G at 

. the present epoch for scales from sub-kpc to ^ 100 Mpc. 

o : 

^ 1 INTRODUCTION 



, Magnetic fields are ubiquitous in the universe and presumably play an important role in most objects in the universe. Their 
' origin however is not well understood (see e.g. Parker 1979; Zeldovich, Ruzmaikin & SokolofT 1983). Galactic fields of micro- 
Gauss strength could have arisen from dynamo amplification of seed fields ~ 10"^" G (see e.g. Ruzmaikin, Shukurov & Sokoloff' 
i-G 1988; Beck et al 1996; Shukurov 2004; Brandenburg & Subramanian 2004). On the other hand, large scale magnetic fields 
could have originated from primordial magnetic fields ~ 10~^ G generated during inflationary epoch in the early universe 
Q ' (Turner & Widrow 1988, Ratra 1992, see Grasso & Rubenstein 2001; Giovannini 2004 for reviews). 

^ , Small seed fields could have arisen from various astrophysical processes during the later stages of evolution of the universe 

(2 < 20) (e.g. Subramanian, Narasimha and Chitre 1995; Kulsrud et al 1997; Grasso & Rubenstein 2001 and Widrow 2002 

. . , for reviews). Harrison (1970) considered a scenario in which a small seed field ~ 10~^^ G is generated owing to the vorticity 

^ ■ in the photon-baryon fluid in the pre-recombination plasma. More recently, Hogan (2000) and Berezhiani & Dolgov (2003) 
^> also considered the generation of magnetic fields from photon pressure in the pre-recombination epoch. From astrophysical 

5_j processes, during the post- recombination universe, the typical scales of the seed field are ~ a fewMpc (see e.g. Subramanian, 
Narasimha and Chitre 1995; Kulsrud et al 1997). On the other hand, magnetic fields can be generated on much larger scales 
in the pre-recombination universe. 

In the pre-recombination universe, photons, electrons, and protons can be treated as tightly coupled fluids. Photons 
however preferentially exert pressure on electrons (the pressure on protons is suppressed by a factor (me/rrip)^). During the 
evolution of the photon-baryon plasma a difference in velocity fields of electrons and protons can therefore be generated, and 
this holds the promise of generating magnetic fields from this induced current. In addition the approximation that photons 
and baryons are tightly coupled, and therefore can be treated as one fluid, breaks down for scales up to several Mpc by the 
time of recombination (during the process of recombination this scale exceeds the horizon scale) . At smaller scales the photons 
free-stream and this in principle can lead to additional contribution to induced currents that might generate magnetic flelds. 

In this paper we study the coupled electron, proton, and photon plasma in first and second order in perturbation theory 
to understand the generation of magnetic fields during the evolution of the plasma. 

In the next section we describe the relevant equations for our study. In §3 we discuss the evolution of magnetic fields 
and its sources in the first and second order in perturbation theory. In §4 discuss our results and give concluding remarks. 
Throughout the paper, numerical values of different quantities are given for the spatially fiat FRW model with = 0.3 and 
r^A = 0.7 (Spergel et al. 2003, Reiss et al. 2004, Tonry et al. 2003, Perlmutter et al. 1999, Riess et al. 1998) with Q,ih^ = 0.02 
(Spergel et al. 2003, Tytler et al. 2000) and h = 0.7 (Freedman et al. 2001). 



2 PRE-RECOMBINATION PLASMA 

The primary components of the plasma in the radiation era which lasts between neutrino decoupling and recombination are 
photons, free electrons and protons. Observations of the Cosmic Microwave Background Radiation (CMBR) which is the relic 
of the radiation existing in this era show that the plasma is almost homogeneous and in thermal equilibrium at the time 
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of recombination (see e.g. Peebles 1993). However, anisotropies observed in the CMBR also indicate that there are spatial 
fluctuations superimposed on this uniform background density. The initial condition for each mode of fluctuation is set before 
it enters the horizon. In the present analysis wo assume that the fluctuations arc adiabatic which means that the entropy per 
fluid particle is conserved. Recent WMAP observations favour this initial condition (Peiris et al. 2003). The electrons interact 
with each other and with protons through Coulomb scattering. The mean free paths for e-e, e-p and p-p collisions are the 
same in this thermal plasma (sec e.g. Shu 1992) and arc much smaller than the astrophysically relevant scales (~ IMpc). 
Hence a continuum description treating them as fluids can be used. In such a macroscopic description the effect of scattering 
between different species is taken into account by including a momentum exchange term in the Euler equation. For photons 
however the dominant interaction is Thomson scattering off free electrons with mean free path (comoving), l^^, at 2 ~ 1000 
for a fully ionized universe being, Z-^e = l/(aaTne) ^ 3Mpc. Here, Ut is the Thomson cross section for e-7 scattering, is 
the electron number density and a is the scale factor. This is comparable to the length scales in consideration and hence a 
Boltzmann particle description is essential. 



2.1 Description of photons 

The photons are described by the phase space distribution function /(x, r;, p, n) where p is the magnitude of photon momentum, 
n is the propagation direction and 77 is the conformal time. Since the distribution is blackbody to zeroth order, we can expand 
it as / = /° + 5f where is the Planck function and 5f is the perturbation. It is convenient to describe the evolution of the 
perturbed distribution in terms of the brightness A(x, r;, n) defined as: 

^=1^ (1) 

The space-time metric in comoving coordinates a;' and conformal time r? for the conformal-Newton gauge is given as (see e.g. 
Ma & Bertschinger 1995 and discussion therein): 

ds"^ = a{r))[-{l + 2^)drf + (1 - 2^)dxidx'] (2) 

Here, a(r;) is the scale factor and 'I', $ are the two potentials characterising scalar perturbations in this gauge. The evolution 
of photons in this metric is then given by the Boltzmann equation for the brightness: 

A + Ti.ftA + n,ft* + $ = C[f] (3) 

Here, over-dots denote derivative with respect to conformal time 77. C[/] is the collision term accounting for the scattering of 
photons with electrons. The linearised collision term for Thomson scattering (neglecting polarisation) is given as (see e.g. Hu 
& White 1997): 

C[f] = UeCTT (^Ao - A + 4ve • n + ^mrijUij'^ (4) 

Here, Ao = J denotes the isotropic part of A and liij is the photon anisotropic stress tensor which takes into account 
the angular dependence of Thomson scattering . It is given by: Hij = J ^ {runj — \5ij) A. The nature of individual terms 
in C[/] suggest that multipole moment expansion in terms of spherical harmonics can give a useful description. The details 
of such an expansion are given in the appendix. We finally arrive at a hierarchy of equations for moment Af^. We can see 
that for a given m, each /-moment A^m is coupled to an Z -I- 1 and an Z — 1 moment in the hierarchy. On the other hand, 
there is no coupling between different m modes. This implies that if there is no source Stm (Appendix A) for a given m and 
if initial conditions are such that A^^ = for that m, then, A^^ = at all times even if other m moments evolve. In the 
present analysis, our emphasis will be on studying the effect of scalar perturbations which correspond to m = 0. For such 
perturbations Hij can be greatly simplified by using azimuthal symmetry about the axis of electron velocity (see e.g. Dodelson 
& Jubas 1994). However our aim here is to study the generation of magnetic field from the evolution of coupled photon-baryon 
plasma. Since the generation of this field can explicitly break this symmetry, one should consider a more general expression 
for the anisotropic stress tensor. 



2.2 Fluid equations for electrons and protons 

As discussed earlier, since the mean free paths of electrons and protons are very small compared to astrophysical scales, we 
can describe their evolution accurately using continuity and Euler equations for an ideal fluid. In linear theory, the density 
fleld Pe,p(x, 77) is expanded as pe,p(x, r}) = pe,p(ri){l + (5e,p(x, rf)), where p is the unperturbed background density and 5 is the 
fractional perturbation. In what follows quantities denoted with a bar on top are background unperturbed quantities. The 
continuity equations for each of the above species is given as (e.g. Ma & Bertschinger 1995): 



(je.p + V • Ve,p - 3i = 



(5) 
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The corresponding Euler equations are: 

d VPe „T aeE ae „ /v-, — Ve\ „ Vj, — Ve 
Ve + -Ve = - V* V, X B + R + (6) 

a pe me me \ T-,e J Tep 

a VPj, _ ^ eE ae „ — Vj, / \ 
Vp + -Vp = ^ - V* + + Vp X B + — ^ — (7) 

a pp nip TTip Tep \mp J 

Here, r^p is the (co- moving) electron-proton coUision time scale; r-ye = l/(neO"Ta) is the photon-electron Thompson scattering 
time scale. E, B are the physical electric and magnetic fields, R = Ap^/Spe, and Pe,p axe the pressures which, for adiabatic 
fluids, can be written as P = P{p)- By taking curl of the Euler equations and using Maxwell's equations, we can get the 
evolution equation for the vorticities {fle.p = V x v^.p) of the fluids as: 



ne + -ne = — — — (aB) Vx(v, xB)+ ^^ ^ 7? 

a mea arj me \ r^e / 

np + -np = — a^B + — Vx VpxB + - — (8) 

a mpU dr) mp AivneeTep \ mp ' 



3 EVOLUTION EQUATION FOR THE MAGNETIC FIELD 

To arrive at the equation governing the evolution of magnetic held, we use the Euler equations for the charged fluids and 
Maxwell's equations (Appendix B). Subtracting Eq. from Eq. (|7J and using Maxwell's equations we first obtain the 
evolution of the current J : 

m„ 9 / J \ d me 1 „ , „. / — v„ \ „m„ m„ J , . 

-f — + = VPc + aE + a Ve x B - R— ^ (9) 

drj XUeJ a e^Ue Uee \ ^le J e UeC'^Tep 

In the above equation we have neglected forces on the proton fluid due to pressure gradient and electric field since they are 
smaller than that for the electron fiuid by the factor me/mp. Taking curl of equation (9) and using Maxwell's equations, we 
get the equation for the generation of magnetic fields: 

1 d , 2■.^^ me „ 9 / J \ me „ / VPe \ „ , me „ / J \ me „ / P(v^ — Ve) \ d me „ / J \ „N 
-j-{aB.) = ^Vx— — +— Vx -Vx VeXB +^Vx U^Vx ^ JLVx — 10 

a ori arj \ne J e \ pe J \neTep J e \ T^,e J a e-' XUeJ 

For studying the generation of magnetic fields in the early universe most of the terms in the above equation can be dropped. 
The first term on the right hand side can be shown to be negligible as compared to the term on the left hand side (see e.g. 
Widrow 2002). Similarly all the terms proportional to B can be neglected if one wishes to study the generation of magnetic 
flelds from zero magnetic fleld initial conditions. These terms can back-react once the magnetic field is generated. We show 
later that the back-reaction terms are negligible for the magnitude of the generated magnetic fleld. These considerations 
simplify the above equation to: 

i|^(a^B) = S(x,^) (11) 
with 

S(x,,) = I^V X f ^\ + ^V X f Mlx^^ (12) 



We now discuss the nature of these source terms of magnetic field generation in first and second order in perturbation theory. 



3.1 Evaluation of the source term: Linear theory 

The source term for any Fourier mode S(k,7)) can be simplified for the linear case. In this case, R = 4p^/(3pe), and r^e = 
l/(ne(jT), are unperturbed quantities and hence don't carry any spatial dependence. The first term of the right hand side of 
Eq. (I12t identically vanishes in this case. The source term can then be written as: 

S(k,7?) = ^^(fl^ -fie) (13) 

Thus, we see that the source for the magnetic fleld in the plasma is the differential vorticity between electrons and photons. 
The vorticity equation for photons is essentially the Boltzmann moment equation for / = l,m = 1 (Eq. IIASfl 'l. We note that 
the source of photon vorticity is 4:Vn/Tje oc fie. This implies that the only source which can excite any Z-moment for m = 1 
is the electron fiuid vorticity. The evolution equation for the electron fiuid vorticity (Eq. Q) shows that the only sources of 
vorticity are the magnetic fleld and the photon vorticity. This means that if the vorticities were zero in the initial condition. 



4 Rajesh Gopal and Shiv K. Sethi 



as is the case with initial zero-vorticity conditions we consider here, none of these quantities can be generated for any scale in 
the linear regime. In particular we can conclude that no magnetic field is generated in linear order for scalar perturbations. 
It should be noted that using Eq. IAS! and Eq. IIA9t allows us to follow modes at which photons are free-streaming at any 
given epoch. Therefore the above conclusion holds for all scales larger than the scales at which electrons and protons can be 
treated as fluids. 



3.2 Source term in the second order 

There are various terms which have to be included in going to second order in perturbation theory. Second order terms can 
arise from treating metric perturbations to second order (Martinez-Gonzalez, Sanz & Silk 1992) or by including the second 
order terms in the electron-photon scattering (Vishniac 1987, Jubas & Dodelson 1995, Hu, Scott, & Silk 1994). Vishniac (1987) 
adopted the simple procedure of including the spatial dependence of densities to include the second order effects. Detailed 
analyses (Jubas & Dodelson 1995, Hu, Scott, & Silk 1994) showed that Vishniac's procedure gives the most important second 
order effect in the electron-proton scattering for sub-horizon scales. This allows us to study scales smaller than the horizon at 
the last scattering surface, ~ 100 Mpc. At larger scales other second order effects from electron-photon scattering and 

the second order metric perturbations might be comparable or dominate. We adopt Vishiniac's procedure here and obtain the 
second order term from treating the spatial dependence of densities i.e. in R, r^e and pe in the source term for magnetic field 
generation (Eq. 112ll 'l. To get estimates of the generated magnetic field we solve for the difference in photon and baryonic 
bulk velocity in the tight-coupling approximation. We argue below that the main contribution to the source ^(k, rj) for any 
scale comes from epochs at which the tight-coupling approximation is valid. The source term is evaluated in the tight-coupling 
approximation in Appendix A and given by Eq. 1A14II . We argue there that for adiabatic evolution, the only term that can 
source the magnetic field generation is given by Eq. IIA15II . This allows us to solve the evolution of the generated magnetic 
field at any scale: 

a'Bik,r,)c^^ J\^'^ J d^fc^(k'),r,')k'xv.((k-k'),V) (14) 

Note that aR is independent of time. Eq. 1 A 1511 can be solved using linear theory evolution of density and velocity perturbations 
for each scale from initial conditions at the time at which the scale is super-horizon to the epoch of recombination, r^rcc 
(Eq. (IA2H '). At the epoch of recombination the photons decouple from the baryons and therefore the source for magnetic field 
generation vanishes. We do not attempt an explicit solution here but seek an approximate understanding of the generated 
magnetic field. We first justify our use of the tight-coupling approximation. As discussed in Appendix A, for each scale the 
tight-coupling approximation is valid for epochs before the Silk damping regime (Eq. iKwi ). From Eq. the magnetic 

field at a given scale k gets contribution from density and velocity perturbations at all scales. It should however be noted 
that if k corresponds to a scale at which the density and velocity perturbations are in damping regime the source of magnetic 
field is also in the damping regime. (More precisely one is interested in the power spectrum of the magnetic field, which, 
from Eq. 11411 . is a four-point function containing density fields. For magnetic field at scale k, the integrand of the source 
is cx P(fc')P(|k — k'l), here P{k) is the power spectrum of the density field; |k — k'| ~ fc if both k and k' are not in the 
damping or free-streaming regime.) This means that most of the contribution to the magnetic field comes from epochs at 
which the tight-coupling approximation is valid. More precisely much of the contribution to magnetic field at a scale k comes 
from epochs rj ^ -qd — i^^^ (Eq. <A20ll '). For scales that are not in the damping regime at -q^^c the upper limit of the integral 
in Eq. 1141 is rj^ac- For smaller scales, the upper limit is ~ i]d- Having identified some generic features of the magnetic field 
source terms we can give an order-of- magnitude estimate of the generated field from Eq. lA15t . For all scales a reasonable 
upper limit on the generated magnetic field for the currently-favoured ACDM model, at a scale L ~ k"^ , is: 

{aB )(L,,7o) = {aB )(L,r,0 ^10 G j (^ lokmsec^i ) [^-) ^^'^ 

Here rji = r^rcc for scales that are not in Silk damping regime at the epoch of last scattering and rji ~ rjd for smaller scales. In 
Eq. 1151 we have used the fact that once the sources of generating magnetic fields vanish, the magnetic field evolves such that 
a^B remains constant (see e.g. Wasserman 1978). In Eq. 1151 . Ve{L) ~ {kP{k)Y^^ (the matter power spectrum is given e.g. 
by Bardeen et al. 1986) and (5e(?7i) is the RMS of the density field. One can compare the magnetic fields at different scales 
by evaluating the sources at the epoch of recombination njiac- As seen from Eq. 1A211 the density field is either non-evolving 
or in the oscillatory phase for much of the period prior to recombination (except for modes r]~^ ^ k ^ (r;/\/3) in the 
matter-dominated epoch). The velocity field however grows oc i] for super-horizon scales in the radiation dominated era and 
for k ^ {t]/V3)~^ in the matter dominated era. Therefore if Eq. I|15|l is evaluated at 77rcc, the small scale magnetic field is 
approximately smaller by a factor (r;cnt/?7rcc); here r^cnt — k~^ is the epoch of horizon entry of the mode k. Using this to write 
the sources of the magnetic field at j^rcc in Eq. 1151 . it is seen that, for all scales, B{r)o) ^ 10"^" G. 
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4 CONCLUSION AND DISCUSSION 

We have studied the possibihty of generating magnetic fields during the evolution of the photon-baryon plasma in the pre- 
recombination universe. For scalar perturbation in linear theory magnetic field is not generated at any scale; this includes 
scales at which the photon-baryon coupling approximation breaks down. We show that in the second order in perturbation 
theory a small magnetic field is generated. The strength of the generated magnetic field is ^ 10~^° for scales from ~ 100 Mpc 
to sub-kpc at the present epoch. 

In Eq. llini, we have neglected several terms which could back-react on the generated magnetic field. It can be readily 
seen that, for the strength of the generated field, these terms are always much smaller than the source term of the magnetic 
field. And therefore we were justified in neglecting those terms for studying the generation of magnetic field. As discussed 
above magnetic fields at small scales are frozen in the plasma from epochs ~ i^J^. It could be asked whether the radiative 
viscosity prior to the recombination can damp these fields. The maximum length scale damped by pre-recombination radiative 
viscosity is oc _B (Jedamzik, Katalinic, & Olinto 1998, Subramanian & Barrow 1998). For the small magnetic fields we obtain, 
the maximum scale of dissipation can be shown to be much smaller than any relevant length scales. 



APPENDIX A: PHOTON BOLTZMAN EQUATION AND INITIAL CONDITIONS 

In this appendix we discuss the complete spherical harmonic moment expansion of the Boltzmann equation without assuming 
azimuthal symmetry of A. The assumption of azimuthal symmetry precludes the existence of vortical modes for the velocity 
fields of photons and electrons. Hence if we are looking for the generation of magnetic fields in the linear order, we have to 
relax this assumption since a possible presence of magnetic fields violates azimuthal symmetry. The notations used in this 
section are self contained. 

The Fourier transformed Boltzmann equation can be written as: 

A + ifc^iA + ifc^i* + $ = — (^Ao - A -I- 4ve • n 4- "^runjllij^ (Al) 
The photon brightness function A can be expanded in terms of the scalar spherical harmonics Yim as: 



A(n) = ^^^i^A,„F^„(n) (A2) 
where, the coefficients Afm are given by the inverse relation. 



= \/^^ / dnYl^{n)A{h) (A3) 



47r 

The photon fiuid variables like over-density 5^ and velocity v-^ are then given by: 

5^= J^A = Aoo (A4) 
. V, = ^ (A5) 

n^ = \kx v^l = ^ (A6) 

(A7) 

By substituting the expansion for A in Eq. lAH and using the familiar properties of spherical harmonics we arrive at the 
following hierarchy of equations for the evolution of the moments Aim- For details of such an expansion we refer to the paper 
(Hu and White 1997). 

A«„ + ik /l^'"" Ae-i^^ + ifc-^^ii^A^+i,™ + = Sim (A8) 

(21 — 1 j (Zl + 6) T^e 



Here, — \/ P — vn?. The source Sim is given as: 

^ 4$ ) 5«5mO + ( ^ - A;*<5,„o ) 5n + A^™5« (A9) 

T-ye ) ) 10 

In the above equations, the coefficients vim are the coefficients in the multipole expansion of the Ve • n term such that 
Ve • n = ^^vtmYtm&n- The vorticity of the photon fluid is tracked by the evolution of the ^ = l,m = 1 moment An. We 
notice from Eq 1A9II . that the source S\\ — Av\i/T.ye This is the only source for the evolution of the m = 1 moment. The first 
two moment equations in the hierarchy give the familiar continuity and Euler equations for photons: 
ii 

(5^ yk ■ Vt - 4$ = (AlO) 
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+ ik^ + n - ik* = ^ (All) 

4 T-,e 

In the above equation, Hi = jikjliij. 



Al Initial Conditions: Tight-coupling approximation 

The electron-proton plasma recombines at a redshift Zrec — 10"^ (see e.g. Peebles 1993). At any epoch in the universe before 
recombination, 77 ^ r^rcc — 2Hq^ (1 + z^cc)~^^'^ /Q^rl"^ , there are roughly five physically relevant length scales: (a) Super-horizon 
scale, k rj~^ (b) scales that are sub-horizon but larger than the sound Horizon, rj/^/i !t k Si r]~^ . At these scales the evolution 
of velocity fields is determined by gravitational potentials, (c) scales smaller than the sound horizon scale but larger than 
the Silk damping scale, rj/y/Z ^ k ^ fcsiik- At these scales the baryon velocity evolution is determined by both gravitational 
potentials and the photon pressure, (d) scales that are in the damping regime but larger than the scales at which photon 
free-stream, fesnk ^ k ^ fefe, fcfs ~ (2Mpc)~^(10^/(l-|-z))~'^. The densities and velocities of baryons decay exponentially in this 
regime (see e.g. Peebles 1980) and (e) k ^ fcts, at these scales photons are free-streaming and therefore photons and baryons 
cannot be treated as coupled fluids. During the evolution in the expanding universe before recombination, the electron velocity 
and density perturbations at most scales first pass through the Silk damping regime before reaching this phase. Therefore 
during this phase 5e,Ve — 0. The only exception to this occurs around the epoch of recombination when the free-streaming 
length increases very rapidly. As the sources of magnetic field generation are nearly zero in this regime, the dynamics of 
plasma at these scales play an unimportant role for our study. In the evolution in linear theory all scales undergo either some 
or all of these phases of evolution. 

Initial condition for each mode is set outside the horizon. Up to phase (c) discussed above, fe <C fefg. During this phase 
the photons are tightly coupled to the baryons and this greatly simplifies the problem (Peebles & Yu 1970, Peebles 1980, 
Hu & Sugiyama 1995). In this approximation, to zeroth order in t^^'- v^, = v^; and 11'-' — 0. Also to this order in r^e'- 
Se = 3/45-,,; this can be readily obtained by subtracting the electron continuity equation from the photon continuity equation 
(Eqs. l|A10|l and lO). To solve for the difference between electron and photon bulk velocity we need to expand to the first 
order in r-,^. To this order, from Eq. l|XTT|l (Peebles & Yu 1970, Peebles 1980): 

1 \ 

^ + -V5-, + Vn + V*j (A12) 

Here the quantities in the bracket on the right hand side are to be evaluated to the zeroth order in Tje- Eq. (IA12II along with 
the evolution of electron velocity (Eq. Q) and the difi'erence of electron and proton velocities (Eq. JUJ) can be used to give 
the following expression for the electric field in the tight coupling approximation: 

aE(x, t)^^( -R^, _ Z£i _ i^V^^^ (A13) 
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In deriving Eq. IIA13t all terms proportional to the magnetic field were dropped as they cannot act as sources for generating 
magnetic field. Taking the curl of this equation and using Maxwell's equation (Eq. iB3l ) one obtains the equation for magnetic 
field generation: 

1 ^.(a^B) = X f-^Ev. + ^ + UvS,] (A14) 



a dr] e \ ^ Pe 

This equation verifies the discussion above that the source of magnetic field generation is electron vorticity in the linear 
theory. With non- vortical initial conditions, Eq. <A14ll shows that all the sources of magnetic field generation are zero in the 
linear perturbation theory. We wish to consider the second order effect by considering the spatial dependence of densities. 
This gives: R = R{Se — 5.y) ~ —l/3RSe in the tight coupling approximation, as 5^ = 3/4(5-y during adiabatic expansion (see 
e.g. Peebles 1980). In second order, the second term on the right hand side is the Biermann battery term. In the adiabatic 
initial condition we consider here, p — p'' with 7 — 5/3. And the source term cx V5e x V5e ~ and therefore in this limit the 
Biermann battery term doesn't contribute. It should be noted that the plasma evolves adiabatically only for scales that are 
not affected by Silk damping (see below). However as the densities and velocities damp in this regime one doesn't expect much 
contribution from these scales. Biermann battery term can also contribute for initial conditions different from the adiabatic 
initial conditions. The third term on the right hand side also vanishes even in the second order in the tightly-coupled regime. 
Therefore the only source of magnetic field generation is the first term on the right hand side of Eq. I|A14^ . Eq. I|A14^ can 
therefore be simplified to: 

-75- a B = — X -5eVe A15 

This equation can be used to get an order-of-magnitude estimate of the generated magnetic field. 

Eq. 1IXT2I can be used to calculate, to the zeroth order, the evolution equation of electron velocity field (Peebles & Yu 
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1970). This equation along with the continuity equation (Eq. @) and V.E — 0, Eq. 1B6L and dropping all terms proportional 
to B, gives: 

■i - ^ k'P^ R k'S-, + ^-^+3<t- (A16) 



a (1 + 7?) p41 + R) 4(1 + -R) ^ a{l + R) 

This equation can be solved along with the evolution equation of S-y by WKB approximation and these solutions can be 
matched to large scale solutions (Hu & Sugiyama 1995). We discuss here approximate solutions at different epochs. First we 
discuss solutions during phase (c) of the evolution. We note that all the terms on the right hand side except for the 5^ term 
are smaller as compared to this term for scales smaller than the sound horizon scale ~ 1/ V^rj. The electron pressure is always 
negligible as compared to the photon pressure in the pre-recombination universe. With these simplification and bearing in 
mind that 1/R <C 1 during the evolution, Eq. I|A16|I is solved to give: 

<S,(k,77) = A(k)cos ( / ujodri'] (A17) 







Here we have only retained the solution compatible with adiabatic initial conditions (see e.g. Hu & Sugiyama 1995) and 
k 



(A18) 



In phase (d) of the evolution of the plasma, the tight-coupling approximation breaks down and the photon-slip which damps 
perturbations (Silk damping) must be taken into account. The solution including the Silk damping is (see e.g. Peebles 1980): 

(5s(k, 77) = ^(k) cos{ujori) exp{—LJdV) (A-19) 

Here, 

^ (A20) 

The silk damping scale, at any epoch, can be obtained from this expression: fcsiik — {15/{2T~feri)y^'^ — (4Mpc)~^((l + 
2:)/10'^)~^''* in the matter-dominated era. The velocity field in the linear evolution remains non-vortical, and hence can be 
found from the continuity equation (Eq. 0). It should be noted that solutions for the baryon density and velocity fields 
differ from the corresponding quantities for the electrons only by replacing R defined here as R' = 4p^/{3pb) (see e.g. Hu 
& Sugiyama 1995). As 7?' 2> 1 for the evolution of the plasma in the pre-recombination universe, for baryonic densities 
compatible with primordial nucleosynthesis, the baryon and electron quantities can be used interchangeably in evaluating the 
second order expression above. 

The evolution of electron density and velocity in the oscillatory regime and the super-horizon solutions, prior to the epoch 
of recombination, can be summarized as (for solutions at super-horizon scales in this conformal-Newton gauge see e.g. Ma & 
Bertschinger 1995): 

(5e oc = constant for k < r]~^ (RD and MD) 

5e = oscillatory for k > ^"^^ (R-D and MD) and for k > rj~^ (RD) 



for rj-^ <k< (MD) 



Ve oc k^fri for k < (RD and MD) and for (^^^ j ~ ^ ~ (MD) 

Ve = oscillatory for k> (^-^^ (RD and MD) (A21) 
Here RD and MD correspond to radiation and matter dominated epochs, respectively. 

APPENDIX B: MAXWELL'S EQUATIONS FOR FRW BACKGROUND 

The Maxwell's equations for the FRW metric in terms of physical fields, E, B, J are as follows: 

V X (a'B) = 4nd'j + (Bl) 
V • B = (B2) 
V X (.^E) ^ (B3) 
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(B4) 



The current J is written in terms of fluid quantities as: 



J = e(npVp — neVe) 



(B5) 



Here, n^.p are the electronic and protonic number densities which are assumed to be equal to the lowest order i.e fie = fip = n. 
From Eq. IBH it follows that V.J = if the second term can be neglected, which is the case here (see e.g. Parker 1979). 
Eq. along with Eq. 1B5II then shows that: 



in the linear theory. 
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